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Abstract 

We study the radiative lepton flavor violating l\ — > Z27 decays in the two Higgs doublet 
model with localized new Higgs doublet in the Randall Sundrum background. We estimate 
the contributions of the KK modes of new Higgs bosons and left (right) handed charged 
lepton doublets (singlets) on the branching ratios of the decays considered. We observe 
that there is an enhancement in the branching ratios with the addition of new Higgs boson 
and lepton KK modes. 



*E-mail address: eiltan@newton.physics.metu.edu.tr 



1 Introduction 



The processes with flavor violation (FV) are worthwhile to study since they exist at least in 
the one loop level in the standard model (SM) and, therefore, they are rich from the theoretical 
point of view. The lepton flavor violating (LFV) interactions are among the most exciting 
candidates of these processes, since they are clean in the sense that they are free from strong 
interactions. Furthermore, the small numerical values of branching ratios (BRs) of LFV decays 
stimulate one to search beyond and to study the more fundamental models in order to enhance 
these numerical values to reach the current experimental upper limits. Among the LFV decays 
the radiative l\ — > ^7 processes reach great interest and their current experimental upper 
limits of the BRs are: BR (// &y) = 1.2 x KT 11 pQ, BR (r -> ej) = 3.9 x KT 7 [2\ and BR 
( r ^/i 7 ) = l.lxl(r 6 (9.0xl(r 8 ;6.8xl(r 8 , 90%CL) [3] (@|; [5] ), respectively. Furthermore, in 
order to search the fi — > e^f decay, a new experiment, aiming to reach a sensitivity of BR~ 10~ 14 , 
at PSI has been described [6J. At present, this experiment (PSI-R-99-05 Experiment) is still 
running in the MEG [7J. 

The theoretical values of the BRs of the radiative LFV decays in the framework of the SM 
are negligible compared to the experimental upper limits and the addition of one more Higgs 
doublet, which drives the flavor changing neutral currents (FCNCs) and the LFV interactions 
at tree level, may cause to pull the theoretical values of the BRs near to the experimental 
upper limits. This is the case that the lepton FV is induced by the internal new neutral Higgs 
bosons, h° and A , and the strength of this violation is regulated by the Yukawa couplings, 
appearing as free parameters which should be restricted by using the experimental data. These 
decays were examined in the framework of the SM with one more Higgs doublet [8]-|llj. the 
so called two Higgs doublet model (2HDM), in [121 H31 HH HSl HS1 [17] . Besides the theoretical 
calculations based on the 2HDM, they were studied in the supersymmetric models [T5]-|24] . 
in a model independent way [25], in the framework of 2HDM and the supersymmetric model 
[26J and, recently, in the SM including effective operators coming from the possible unparticle 
effects [27]- [28]. 

Another possibility to enhance the numerical values of the BRs of these processes is to 
consider the extra dimension which results in the additional effects of the KK modes of the 
particles in the loops, after the compactification. In the present work, we consider the extended 
Higgs sector, the 2HDM, in the the RSI background [2HJ ED]- The RSI model is based on the 
curved extra dimension and the corresponding metric reads 

ds 2 = e' 2a r ] ^dx tl dx u - dy 2 , (1) 
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where a = k\y\, k is the bulk curvature constant, the exponential e~ CT , with y = R\6\, is 
the warp factor. Here, the extra dimension is compactified onto S 1 jZi orbifold and R is the 
compactification radius. The extra dimension has two boundaries, the hidden (Planck) brane 
and the visible (TeV) brane, with opposite and equal tensions. This choice leads to the fact that 
the low energy effective theory has flat 4D spacetime, even if the 5D cosmological constant is non 
vanishing. In the RSI background, the gravity is taken to be localized on the hidden brane and 
to be extended into the bulk with varying strength and the SM fields live in the visible brane. If 
some of the SM fields are accessible to the extra dimension, the phenomenology becomes richer 
and there are various work done in the literature respecting such scenarios [31]- [H]. If fermions 
are accessible to the extra dimension and there is a Dirac mass term in the Lagrangian, the 
fermion mass hierarchy can be explained. In this case the fermion mass hierarchy comes from 
the possible fermion field locations [35j EH USE HQ] . The quark and lepton FV, which is based 
on the different locations of the fermion fields in the extra dimension, is extensively studied in 
[45. 46J. In these works, it is considered that the FV is carried by the Yukawa interactions, 
coming from the SM Higgs-fermion-fermion vertices. In [47], the high precision measurements 
of top pair production at the ILC is addressed by considering that the fermions are localized 
in the bulk of RSI background. In recent works j48j H9], the various experimental FCNC 
constraints and the electro weak precision tests for the location parameters of the fermions in 
the extra dimension are discussed. The other possibility is to consider the localization of Higgs 
field in the extra dimension. The brane localized mass terms for scalar fields are considered in 
order to get small couplings of KK modes with the boundaries [35] and these mass terms result 
in that the zero mode localized solution is obtained. In [41] the hierarchy of fermion masses is 
analyzed by taking that the Higgs field has an exponential profile around the TeV brane. [42J 
is devoted to an extensive work on the bulk fields in various multi-brane models. 

In our work, we assume that the new Higgs doublet is accessible to the extra dimension 
of RSI background. First, we study the case that the charged leptons are restricted to the 
4D brane and, second, we consider that the charged leptons are also localized in the extra 
dimension. Notice that, in both cases, the gauge bosons are necessarily accessible to the extra 
dimension. 

The paper is organized as follows: In Section 2, we present the BRs of LFV interactions 
h — > hi in the 2HDM, by considering that the new Higgs doublet is localized in the extra 
dimension of RSI background. Section 3 is devoted to discussion and our conclusions. In 
Appendix A, we study the construction of new Higgs boson mass matrix. In Appendix B, we 
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present the amplitudes appearing in the calculation of the decay widths of the radiative decays 
under consideration Appendix C is devoted to calculation of the zero mode lepton fields and 
their KK modes. 



2 LFV li — > Z27 decays in the Randall Sundrum back- 
ground with localized new Higgs boson 

We start with the action for the new Higgs doublet (see for example [31], (12] for a massive 
bulk scalar field case), 

= ~ fd'xfdy^ (g MN (cW)t d N + ml 4>) , (2) 
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where g = Det[gMN] = e~ 8a , M,N = 0, 1, ...,4. The decomposition of the scalar doublet into 
KK modes 

00 

<}>(x,y) = Y,<t> (n) (x)fn(y) (3) 

n=0 

brings the action eq.(T5]) into form 

5s = -E / d*x [rT (dMx)^ d v (f> n (x) + m 2 nS (</> n (x))t n (x)) , (4) 

n=0 

with the second order differential equation 

- e " 17, ( e_4ff dJ i !l ) + m l = <8*° /»(*) , (5) 



and the orthogonality relation 

R 



dye- 2 °f*(y)f m (y) = 5 nm . (6) 

ttR 



The choice of the mass term in eq.© 



,dcr si .da 



results in the differential equation in the bulk 



- e 4CT -f (e" 4 " + a k * fn{y) = m \ s <?° f n (y) , (8) 



dy \ dy 

where n = 1,2, Now, the the boundary mass ter 

cP_a 
dy 



m \,bound — b 1 2 ' 



1 Here the boundary mass terms have the same magnitude and the opposite sign on the branes. The idea of 
brane localized mass terms has been considered for scalar fields in [31j . |35) . 
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is considered in order to obtain zero mode Higgs doublet and this term induces the boundary 
condition 

(^^-bk<P(x,y)y =0 ^ R = 0. (10) 

Notice that the non- vanishing zero mode can be obtained with the fine tuning of the parameter 
b and a, 



b = 2 + V4 + a, (11) 

and it reads 

bky 

f (y) = , . (12) 

' e 2(b-l)kitR_l 



(6-1) k 

On the other hand, the KK mode Higgs doublet is obtained as 

fn(y) = ^— ^v^(e ff x nS ) + a n Y VI ^(e a x nS )^j , (13) 
where Ng n is the normalization constant, x n s = and a n reads 

( 2 - b) J^E(x nS ) + X nS ^/j+^Qns) 

ttn (b-2)Y v ^ E (x nS )-x nS Y^(x nS )- { } 

Here, the functions Jp{w) and Yp(w) are the Bessel function of the first kind and of the second 
kind, respectively. Finally, the mass spectrum of KK modes (n = 1, 2, ...) are obtained by using 
the boundary conditions at y — and y = n R (see eq.( [ID]) ). 

m n5 ~ (n + \ (b - 2) - |) vr fc e~ fc ^ , (15) 

for ke~ knR < m nS < fc. 

At this stage, we consider two possibilities for the charged leptons : 

• they are restricted to the 4D brane 

• they are localized in the extra dimension. 



2 There is another possibility of fine tuning of the parameters b and a for the non-vanishing zero mode, 
namely 6 = 2 — vi + a. However we ignore this choice since it is not appropriate for the brane localized fcrmion 
scenario and bulk fermion scenario with the parameter set used (see discussion section for details). 
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2.1 The charged leptons restricted to the brane 

The LFV interactions are driven by the part of the action 



s Y = / d b x^ lufaEjR + h - c - 6 (y - nR ) 



(16) 



where L and R denote chiral projections L(R) = 1/2(1 =1= 75), 2 is the new scalar doublet, In 
(Ejn) are lepton doublets (singlets), ^f^-, with family indices i,j , are the Yukawa couplings 
in five dimensions, which are responsible for the flavor violating interactions in the leptonic 
sector. Here, we assume that the Higgs doublet 0i lives on the visible brane and it has non- 
zero vacuum expectation value in order to ensure the ordinary masses of the gauge fields and 
the fermions. On the other hand the second doublet, which is accessible to the extra dimension, 
has no vacuum expectation valuqj: 



1 

V2 







+ 



V2X + 
iX° 



and 



< 01 >= 



wo 

y/2 V v 



1 

V2 



< 02 >= 



V2H+ 
H 1 + iH 2 



(17) 



;is) 



The new Higgs doublet 02 is expanded into KK modes after the compactification of the extra 
dimension as given in eq.((3]) and the zero (KK) mode Higgs fields are obtained by imposing the 
fine tuning condition in eq. (fTT|) . On the other hand, after the electro weak breaking, the SM 
Higgs acquires the vacuum expectation value eq. (|T8l) and there appears mixing between zero 
mode and KK mode new Higgs bosons. However, we do not take into account the additional 
effects coming from this mixing since they are suppressed (see Appendix A for detail). 

For the effective Yukawa coupling we integrate out the Yukawa interaction eq. (}T6l) over 
the fifth dimension by taking the zero mode neutral Higgs fields S = h°, A : 



where 



Vu 



wR 



-ttR 



ij S 5 ' 



(19) 



dye-^f (y)5(y-7rR) 



Here we take the Higgs potential as 



Cl 



0! - v 2 /2f + c 2 {4>i^f + c^i^+fo) - (<t>tfa)(<ti<h)] + c 4 [i?e(0+0 2 )] 2 + c 5 [/m(^+0 2 )] 5 



This choice leads to no tree level mixing between the CP even neutral Higgs bosons, namely H° and H . 
Therefore, the SM particles (new particles) are collected in the first (second) doublet and H 1 , H 2 are obtained 
as the mass eigenstates h° and A respectively. Notice that, in general, the mixing between the CP even neutral 
Higgs bosons can exist in the loop level when one considers the quantum corrections. 
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Jb-i) kwR 

• (20) 

e 2 (6-1) k 7T R._i 
fc(b-l) 

Here it is assumed that the coupling • is flavor dependent and it is regulated in such a way 



that the overall quantity ■ is pointed to the chosen numerical value of £f 



>5ij 

'ij^5ij 10 P^inteti uu uiic ^iiuocii iiuiiiciil,cxi vaiuc ui t,^-. 

The effective Yukawa coupling inducing the tree level interaction among the KK mode Higgs 
and charged leptons reads 



t5 n = V3 n & v , (21) 



where 

V t = e -TT— U-2(e k * R x nS ) + a n Y b _ 2 (e k * R x nS )) . (22) 

Here N Sn is the normalization constant (see eq.( lT3|) ) and b (a n ) is defined in eqs.f JTT]) (([1 
Finally, the effective Yukawa coupling is obtained as 



Now, we present the decay widths of the LFV Zi — >• Z27 decays, including the KK modes of 
new neutral Higgs fields. Since these decays exist at least in the one loop level, there appear 
the logarithmic divergences in the calculations. In order to eliminate these divergences, we 
follow the on-shell renormalization scheme. In this scheme, the self energy diagrams can be 
written in the form J2(p) — iP~ m h)^2(.P)(,P~ m h) 5 which results in that these diagrams do not 
contribute for on-shell leptons and, only, the vertex diagrams (see FigJT]) contributcl. Taking 
only tau lepton for the internal linej, the decay width V reads 

T(h^l 2l ) = c 1 (\A 1 \ 2 + \A 2 \ 2 ) , (24) 

where 

Ax = A\ + Af KK , 

A 2 = A° 2 + Al KK . (25) 
For the explicit expression of these amplitudes see Appendix B. 



4 This is the case that the divergences can be eliminated by introducing a counter term with the relation 
yRen _ yo _|_ yC ^ w h ere yRen ^yo^ j g ^ ne renormalized (bare) vertex and by using the gauge invariance 
fcnyRen _ q jj er6; fefi j g fo e f our momentum vector of the outgoing photon. 

5 We take into account only the internal tau lepton contribution since we respect the idea that the couplings 
$Nij (*'■? = C'l 1 )! are small compared to Ti (i = e,/i, r), due to the possible proportionality of them to the 
masses of leptons under consideration in the vertices. Here, we use the dimensionful coupling £jy f • with the 
definition { - = J £n ij where N denotes the word " neutral" . 



2.2 The charged leptons localized in the extra dimension 

The part of the action which drives the LFV interactions in this case reads 



Sy = J d 5 X^ Ugj l iL <P 2 E jR + h.c\ 



(26) 



where In {Ejp) are lepton doublets (singlets) which are localized in the extra dimension. The 
addition of Dirac mass term to the lagrangian of bulk fermions causes this localization [32| IM| 
|35| [37] |38| HOI HI]. Since the combination ipifj is odd due to the two possible transformation 
properties of fermions under the orbifold Z 2 symmetry, Z 2 i\) = ±75'?/', in order to construct the 
Z 2 invariant mass term, one needs Z 2 odd scalar field to be coupled. This discussion leads to 
the mass term 

S m = -j^ji vrgm{ y)W, (27) 

where m(y) = m^-^- with cr'(y) = jj^. With the help of the given mass term the localized 
zero mode leptons are obtained. We present the construction of the zero mode and KK mode 
leptons in the Appendix C extensively. 

For the effective Yukawa coupling similar to the previous case, we integrate out the 
Yukawa interaction eq. (I2"6"j) over the fifth dimension. By taking the zero mode lepton doublets, 
singlets (see eq. (1571) ) and neutral Higgs fields S = h°,A° (see eq. (fT2]) ). we get 



iij ((£y) ) — Vs°RL(LR)ij %5ij ((^5ij) ) ' (^8) 



where 



ttR 



v SRLij = 7T~5 / dy XiRo(y) Xjw(y) fo(y) 

2 7T H J-ttR 

M _ e (b-r iR -r jL )kTr R^j 



(r- J_-r h\ j l-e^- 2r rR)k^R / l- e ^~ 2 r j l) k * R / l_ e 2 (b-1) fcTT fl 

\r iR + r jL 0)\j {2 r iR -l) V (2r jL -l) V k (1-b) 



(29) 



Here, similar to the previous scenario, the coupling in five dimension is flavor dependent 
and it is regulated in such a way that the overall quantity V RL , LR \ ^ £f^- is pointed to the chosen 
numerical value of £y f(£|/)M- This is the case that the hierarchy of new Yukawa couplings, 
describing the tree level Higgs zero mode(5'^ - ) )-lepton zero mode (/^-lepton zero mode (l^) 
interaction, is not related to the Higgs field and lepton field locations. 

The effective Yukawa coupling which drives the tree level KK mode Higgs (S^ n ')-l^-l^ ' 
interaction is 

Zg n ° ((4? n °) f ) = Vf RL{LR)ij ^ ((^)t) , (30) 



where 

yn0 = I-nR dyeV- r **- r ^° (J b -2(e°x ns ) + a n Y b ^x nS ) J 

iV -5ny k(2r iR -l) V *<?r jL -l) 

Using the eqs. ( 1301) and ( 13TT) . the effective Yukawa coupling £^ n0 is obtained as 



VS% r .. (n R + r jL -b)J ^^ 



X 



jT* eC 2 -^"^) CT ( J,„ 2 (e CT x nS ) + a n Y fe _ 2 (e CT x nS )) • (32) 



3(0) — /(°) — /(") vertex drives another possible tree level interaction appearing in the loop 
calculations and the corresponding is effective Yukawa coupling reads 

Sij I ) — V SRL(LR)ij^5ij \\$5ij) J ' 

where 

N Ln{Rn) Z*r dye^wM)" (ji^ w (e , ^) + ci W yi TrjlW) (fi'^) 



l/ 0n 

V SRL (LR)ij 



p / 1 _ e ( 1_2r iii(iL)) fc ' ri? e 2[b-l)kiT R_i 



Here the parameters x n R(L), cr{L), the lepton localization parameters r^/jn and the normal- 
ization constant N R ^ n are given in Appendix C. By using the eqs. (1331) and ( 1341) we get the 
effective Yukawa coupling £^ 0n ((£•■• 0n )"M as 



tEQn((tE0n\\\ _ *SHHLR)ij t E ( ( t B\\\ _ AT (1 ~ e (1 2 r * L Q«>> k * R ) 



(34) 



2 : F''jZ,(jfl) 



[e a XnL(R)) 



(iL) -rjt (jR) )knR 



(35) 



Finally, the tree level - - interaction is carried by the effective Yukawa coupling 
ifj mn and it reads 

tfr n ((e| mn ) f ) = v s T L {lr) y ef« , (36) 



8 



with 



V3& {Lm = NLn Tl m ^ dye^- r ^°(j b - 2 (e°x mS )+a n Y b „ 2 (e°x mS ) 

X ( ^t^m^^^l^)) + c L(i? )Fi Trji(jfl) (e CT a; nL(R) ) ) , 



(37) 



and, therefore, we get 



\\ _ S RL {LR) ij C E f(tE\\\ 



Y T fl(2r ii(ifl) -1) Y CS=05 ^ ( iL ) + ~ °> 

N Sm (l — e ( b - r iR(iL)- r 3 LUR)) k7rR ^ 

x r R dyed-'-iHw)^ f J 6 _ 2 (e CT x mS ) + a n F 6 _ 2 (e (T x m5 ) J 



(38) 



The decay widths of the LFV fi — > ej, r — > and r — > p,j decays are calculated by using 
on-shell renormalization scheme (see the section I2.ip and we have 

r(/ 1 ^/ 27 ) = Cl (|A 1 | 2 + |A 2 | 2 ), (39) 

with 

A 1 = A° + A? KK + A[ KK + A(> lKK , 

A 2 = Al + A s 2 KK + A l 2 KK + Al' lKK . (40) 

Notice that we present the explicit expression of the amplitudes A®^, Af&j? and Af'}Jf K in 
Appendix B. 



9 



3 Discussion 



The Yukawa interactions coming from lepton-lepton-5' vertices drive the radiative LFV l\ — > Z27 
decayqj, and their strengths are regulated by the Yukawa couplings which are free parameters of 
the model used. In the present work, we study these LFV decays in the RSI background and we 
assume that the new Higgs doublet and the gauge fields are accessible to the extra dimension. 
Here, in order to obtain zero mode Higgs doublet, one considers the boundary mass term (see 



eq.®) and impose the fine tuning 6 = 2 + a/4 + a (eq. ffTT]) ) of the parameters 6 and a which 
regulates to the boundary and bulk mass terms. Finally, the zero mode new Higgs doublet 
is obtained as an exponential function of the parameter 6 (eq.( lT2"j) ) and it is highly localized 



around the visible brane. The choice 6 = 2 — \/A + a is also possible for the non-vanishing 
zero mode. However, we do not take this possibility into account because of the following 
reason: If the fermions are localized on the 4D brane, the overall quantity V^f^- is fixed to 
the chosen numerical value of with the assumption that the coupling ^ ■ is flavor dependent 
and appropriately regulated. For 6 = 2 — \/A + a the coupling V$ is a number of orders smaller 



compared to 6 = 2 + a/4 + a and, since this term appear in the denominator of the coupling 
£ En according to our definition, £ n exceeds the range of perturbative calculation. Notice that 
the coupling V$ n is not so much sensitive to the parameter 6. In the case of bulk fermions, the 
choice 6 = 2 — \/A + a results in extremely large coupling which breaks the perturbative upper 
limit (negligible coupling which causes weak sensitivity to KK mode contributions) for set II 
(set I) . 

As a first attempt we assume that the leptons are restricted to the 4D brane. In this 
case the contribution of the extra dimension is due to the new Higgs KK modes which appear 
in the internal line of the loop with the modified Yukawa couplings (eq. (1231) ) . Second, we 
consider that the leptons are also localized in the extra dimension. We follow the idea that 
the localization of the lepton fields in the extra dimension occurs with the addition of a Dirac 
mass term mi = ra' with o = k\y\ (eq. (j2"7|) ) . In this case, the right and left handed lepton zero 
modes (eq. (|57p ) are chosen to locate at different positions in the extra dimension in order to 
explain different flavor mass hierarchy. In the scenario we choose the contribution of the extra 
dimension is coming from the new Higgs KK modes and the lepton KK modes appearing in 
the internal line of the loop. The FV is carried by the new Yukawa couplings which are fixed 



6 Here, we do not take into account the internal neutrino mediation due to their weak contribution to the 
BRs of the processes we study and, therefore, we assume that the lepton FV comes from the internal new 
neutral Higgs bosons, h° and A . Notice that we ignored the possible restrictions coming from the hadronic 
decays. 
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to an appropriate number, respecting the current measurements and the location parameters 
of leptons are responsible for the lepton mass hierarchy. This choice makes the constraints 
coming from various LFV processes to be more relaxed. Here, we consider two different set 
of locations of charged leptons in order to obtain the masses of different flavors^. In the first 
set (Table [2]), we consider the left and right handed fields having the same location in the 
extra dimension. In the second, we choose the left handed charged lepton locations as the 
same for each flavor, and we estimate the right handed ones by respecting the current charged 
lepton masses. For the second set, we observe that the BRs of the decays under consideration 
enhance since the KK mode couplings to the new Higgs scalars, which are highly localized near 
the visible brane, become stronger if the left handed lepton field is near to this brane. For 
the effective Yukawa couplings in four dimension we choose that h3 = e ;/ i are smaller 
compared to £,§ Ti i = e, //, r, since latter ones contain heavy flavor and we assume that, in four 
dimensions, the couplings is symmetric with respect to the indices i and j. Furthermore, 
the curvature parameter k and the compactification radius R are among the free parameters of 
the theory. Here, we take k R = 10.83 and consider in the region 10 17 (GeV) < k < 10 18 (GeV) 
(see the discussion in Appendix C and [i0]). Throughout our calculations we use the input 
values given in Table (OQ). 



Parameter 


Value 


m T 
m h o 
m^a 
Gp 


0.106 (GeV) 

1.78 (GeV) 

100 (GeV) 

200 (GeV) 

1.1663710- 5 (GW- 2 ) 



Table 1: The values of the input parameters used in the numerical calculations. 

In the case that the leptons live on the 4D brane the contribution of the Higgs boson KK 
modes is negligible for the decays under consideration The weakness of the new contribution is 
due to the tiny ratio zs n appearing in the expression (eq.( l47|) ) which represents the additional 
effects to the amplitudes. 

Now we analyze the case that the leptons are also accessible to the extra dimension. 

The gauge sector is necessarily lives in the extra dimension and their KK modes appear after the compact- 
ification of the extra dimension. The different fermion locations can induce additional FCNC effects at tree 
level due to the couplings of neutral gauge KK modes-leptons and they should be suppressed even for low KK 
masses, by choosing the location parameters tl (j~r) appropriately. In the set of location parameters we use 
(Tabled]), we verify the various experimental FCNC constraints with KK neutral gauge boson masses as low as 
few TeVs (see the similar the set of location parameters and the discussion given in [48l 149] .) 



11 



FigfS] represents the parameter k dependence of the BR of the LFV \i — > e^y decay for 
£,N,re = 0.01 GeV, C^ru = 10 GeV. Here the solid (dashed, short dashed) line represents the 
BR without KK modes of leptons and new Higgs bosons (with KK modes of leptons and new 
Higgs bosons for lepton location set II, set I), for a = 0.01 and O.lo It is observed that the 
BR (/i — > e^y) is of the order of 10~ n without the internal lepton and new Higgs boson KK 
mode contributions. The addition of these KK modes result in that the BR enhances almost 
2x one order, for the lepton location set II, especially for the small values of the parameter k. 
For the set I, the enhancement of the BR is negligible. On the other hand, the BRs are weakly 
sensitive to the parameter a which plays a crucial role in the localization of new Higgs bosons. 
We present the parameter a dependence of the BR (/i — > ej) for the lepton location set II, in 
Figj3]for k = 10 18 GeV. This figure shows that the enhancement of the BR (/i — > ey) is of the 
order of ~ 0.1% in the range of a, 0.01 < a < 1.0. This is a negligible enhancement which can 
not be determined. On the other hand the enhancement in the case of set II is due to the fact 
that the left handed leptons (KK modes) are near to the visible brane and their couplings to 
the new Higgs bosons, which are localized near the visible brane, become stronger. 

FigJH is devoted to the parameter k dependence of the BR of the LFV r — > cy decay for 
ifi Te = 0.1 GeV, ifi TT = 50 GeV. Here the solid (dashed, short dashed) line represents the 
BR without KK modes of leptons and new Higgs bosons (with KK modes of leptons and new 
Higgs bosons for lepton location set II, set I), for a = 0.01 and 0.1. This figure shows that the 
BR (r — > &y) is of the order of 10 -12 without the internal lepton and new Higgs boson KK 
mode contributions. The addition of these KK modes results in that the BR enhances almost 
three orders for the small values of the parameter k and the lepton location set II. For the 
set I, the BR enhances to the value almost two times larger compared to the one without KK 
modes. Figj5] represents the parameter a dependence of the BR (r — > ej) for k = 10 17 GeV, for 
the lepton location set II. It is observed that the enhancement of the BR (r — > cy) is greater 
than ~ 2.0% in the range of a, 0.01 < a < 1.0. Similar to the previous decay, this is a small 
enhancement which can not be determined. 

Figj6] shows is the parameter k dependence of the BR of the LFV r -> /ry decay for 
£n,t[i = 1-0 GeV, £Jv,tt = 50 GeV. Here the solid (dashed, short dashed) line represents the 
BR without KK modes of leptons and new Higgs bosons (with KK modes of leptons and new 
Higgs bosons for lepton location set II, set I), for a = 0.01 and 0.1. The BR (r — > fi'y) is 
at the order of the magnitude of 10~ 10 without the internal lepton and new Higgs boson KK 



For a = 0.01 and a = 0.1 the curves almost coincide 
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mode contributions. The lepton and Higgs boson KK modes cause more than three order 
enhancement in the BR for the small values of the parameter k and for the lepton location set 
II. For the set I, this enhancement is almost two times of the BR without KK modes. In FigJT] 
we present the parameter a dependence of the BR (r — > /17) for k = 10 GeV and the lepton 
location set II. We observe that the enhancement of the BR (r — > ^7) is more than ~ 2.0% in 
the range of a, 0.01 < a < 1.0. 

At this stage, we would like to summarize our results. For the brane leptons, the contribution 
of the Higgs boson KK modes to the BRs of the radiative LFV decays is too small to be detected. 
However, if one considers that the leptons are also accessible to the extra dimension, there exists 
a considerable enhancement in the BRs, especially for the small values of the parameter k. This 
enhancement occurs for the lepton location set II and it is due to the fact that the left handed 
leptons (KK modes), which are near to the visible brane, have enhanced couplings to the new 
Higgs bosons, which are also localized near the visible brane. Finally, we observe that, the 
BRs are weakly sensitive to the parameter a which regulates the localization of the new Higgs 
doublet in the extra dimension. With the more accurate forthcoming measurements of the BRs 
of the LFV decays it would be possible to test the existence of the warped extra dimensions 
and, to get a considerable information which fields are accessible to the extra dimension. 
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Appendix 



A The mass matrix of new Higgs boson 

The Higgs potential which creates the masses of neutral CP even and CP odd Higgs bosons, 
S = h°, A reads 



V, 



Ills 



c / ,[i?e(0+0 2 )] 2 + c / A [/m(0+0 2 )] ; 



(41) 



where (f)% and 02 are given in (eq. (|T7|) ). After the electrowek breaking the SM Higgs acquires a 
vacuum expectation value (see (eq. (118jl ) and the mass lagrangian of new CP even Higgs boson 
(h°) § becomes 

2 

I , (42) 



£s = ~ E mis (S^(x)f +c' h V -( S<®{x) + E S™{x) a n ) , 

Z n=\ Z \ n=l / 

with Ch = d h /o (it R), a n = ^^ji and m| = Ch v 2 , S = h°. By using the mass lagrangian the S 
boson mass matrix is obtained as (see [50] and [51] for boson mass matrix in the one and two 



non-universal extra dimensions, [SB] for U(l)y gauge boson mass matrix. 



Ml 



( m 2 s 
«i m| m\ s + m 2 s a\ 



0(2 m 2 s 



«2 «i Tris 



«2 Trig 
m 2 s oi\ «2 
ml s + m 2 s oi\ 



\ 



\ 



(43) 



and the determinant equation reads 



det (Mg -XI 



n { ni 



2 

n S 



(XJ 

A)) (m 2 s - A - A m| E 



a, 



n=l n=l m nS 

This equation is used to calculate the physical masses of zero and KK modes of S bosons and 
their eigenstates. To leading order, the physical S boson mass (the zero mode one) reads 



A 



0. 



(44) 



phys \ 2 



i + E 



m| a 2 n 



2 / ' 

n=l m nS ) 



(45) 



since m n s » ms- Notice that in our numerical calculations we do not take into account the 
additional effects coming from the mixing because they are suppressed due to the fact that the 
KK mode masses are considerably larger compared to the zero mode one. 



3 The similar mass lagrangian appears for the CP odd Higgs boson A with the replacement c' h — > c' A 
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B The amplitudes appearing in the text 



Here, we present the amplitudes which appear in the calculation of the decay widths of the 
radiative decays under consideration. 

In the case that the charged leptons are restricted to the brane, the amplitudes A®, A®, 
Af KK and Af KK in eq.(|25]) read 

< - «^{ m " ^ c r * * <* + » - ^ + ^> 

- ^(LrSir / dx I dyxy{^- + 



o 



+ m r ^N,rh £n,Ht J q dx J Q d V ( X ~ X ) ~ ^)| > 

JO Jo LfrQ L A o 

- ^^^/^^(i-l)^-^, (46) 
^ ^ = QrT^j E K ^ S,?h * (* + V - 1) (f^ + 

4m r n^l I JO JO L h n L A n 

~ m h^N,i 2 rU,hr L dx j dyxy{- + 



'o Jo L h o n L A o 

+ rn T IZi, iNlr Jl dx f*~* dy(x-l)( Zh ° n * A ° n 



n 



Lh° n L A o n J 

Al KK = Qrj^ £ l-m h ~e N l T el? lT f dx f~ X dyx(x + y-l)(f^ + ^) 



E ,, TE n I dx I J ' n ' n 

) Jo 

1 i-l-x 



o Jo L h o n L A o 



n 



- m T i E N l T g» x / dx f dy (x - 1) (f^ - , (47) 

JO Jo L h o n L A a n 



where 



L s = z s + x 2 z s + x(l + (y-2)z s j 
L Sn = z Sn + x 2 zsn + x{l + (y-2)z S n) , (48) 

with z s = ^f-, z Sn = Here , h (l 2 ) = T;/J,(fjL or e;e), c x = Gf ^ZT H , A x (A 2 ) is the left 

S n S 

(right) chiral amplitude, Q T is the charge of tau lepton and m n s is the internal Higgs KK mode 
mass (see eq.(IT5l)). Notice that we take the Yukawa couplings real. 



15 



If the charged leptons are also accessible to the extra dimension the amplitudes A\ , 2 % , Af t 
and A x \ 2 s (see eq. pojl ) are 

Al = QT-^\m h £§, Tl J§, Th f dx f X dy X (x + y-l)(^ + ^) 
4 m; ( ' ' Jo Jo L h o L A o 

~ m i2^NMT^N,hr f dx f dyxy{^- + 

JO JO Liyfi ■L'A 

+ m r ^n,t1 2 €wit J dx J dy(x- 1) (|^- - -^-) j , 
A ° = «-4^{- m -^® A ,/>/"*x(« + v-l)(^ + ^) 
+ m h ^N,ri 2 %N,Th J dx J dyxy(j^-+ ■ |^-) 

- rrir^i ^§ h f dx f dy(x- 1) - ^-)\ , ( 
Jo Jo Lyfi L A o J 



A? KK = Qrj^ t {m h m% [ dx r* dyx(x + y-l)(^ + ^) 

4m r „=i 1 Jo J o L h o n L A o n 

FEnO FEnO f 1 j [ * J / z h° n Z A° n n 

- mi 2 ^ N .i 2 r^N,hr / dx / dy X y (- + - ) 

JO JO L h o n L A o n 

+ rn T t™ f dx t X dy (x - 1) (f^ - , 

Jo jo L h o n L A o n J 

4 M = Qr^~ 2 t\-m h ,t /' dy x (x + y - 1) (f^ + ^) 

™; „tl { Jo Jo L h o n L A o n 

FEnO FEnO f J /" X J / Z h°n . Z A° n \ 

+ m hU,ri 2 U,rh / ^ / dyxy{- + - ) 

JO JO I^o n L A o n 

'Wo Jo L h o n L A o n J 



48 r< ^ { m 2 nR V 

2 ^ 

+ ~ 2 1- m; 2 £w,?Tt (p( z nL,h°) + G( z nL,A< 



48 mf m: - v 



r n=l I "'nil 

2 



A S,IKK Qr \ -> J m r „ / c £ mn\ t f <F-Emn\t /Vv -v ^ I /^/"~ 

48 m r „, m =l l m ni? 
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o 

) + G{z nLy 

m nL ) 

a s,ikk —Qt \ - J m T m i c E mn \\ i c E mn \\ ( n i \ ^n(^ 

48 m r n,m=l { m nR 

+ ^m l J^^(G(z nL , h0m ) + G(z nL , A o m ))\. (52) 



,,4° m 



Here Ai (A 2 ) is the left (right) chiral amplitude, l\ (I2) = r; n {p, or e; e), the functions F(w 
G(w) are 

w (3 — Aw + w 2 + 2lnw) 



F(w) 



[l-w 



. w (2 + 3w — 6w 2 + w 3 + 6wlnw) , . 

G(w) = o^? ' (53) 

Gf.a em mf m 2 . _ mg _ m \h (nfl) _, _ m n£ (nR) 1 r f 

C l - 32^ ' ' S ~rr& ' ^ 72 r7?T' Z nL(nR),S ~ m * , ZnL{nR),Sm ~ m ± WltH Lett 

(right) handed internal lepton KK mode mass m n L( n R) (eq.( |68|) ). In eqs. ( H9|) and (!50|) the 
functions L5 and Ls„ are given in (eq.( )48p ). In the case that the incoming and outgoing lepton 
masses are ignored in the functions L s and L Sn one gets the integrated form of ^4?( 2 ) as 

(54) 

Notice that, for the amplitudes an d AffJ? K , the incoming and outgoing lepton masses 

are ignored in the functions Ls and Ls n since the internal KK leptons are heavy. 

C The construction of zero mode and KK mode leptons 

This Appendix is devoted to the construction of the zero mode and KK mode leptons in the 
case that the leptons are localized in the extra dimension with the help of the Dirac mass term 
given in (eq. ( 12~T|) ) . We start with the expansion of the bulk fermion as 



oc 



4>(x»,y) = -_ £ ^ n \x»)e 2 ° Xn {y). (55) 
V Ztt K n=0 



By using the normalization 

1 r R 



2nR J-ttR 
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dy e a Xn(y) Xm(y) = S nm , (56) 



and the Dirac equation the zero mode fermion is obtained as 



Xo (y) = N e- r % (57) 
where r = m/k and Nq is the normalization constant: 



knR(l-2r) 

^o = y efc ^ (1 _ 2 ,)_r ( 58 ) 

The appropriately normalized solution 

x' {y) = e" f Xo(y) (59) 

is localized in the extra dimension where the localization is regulated by the parameter r. For 
r > | (r < |) this solution is localized near the hidden (visible) brane and it has a constant 
profile for r = ~. 

Now, we are ready construct the SM leptons. What we need is to consider SU{2)l doublet 
ipL and singlet ipR with separate Z 2 projection conditions: Z 2 ipR = ■jn'ipR an d Z 2 ipL = — l^ipL 
(see for example [32J ) . Finally we get the leptons accessible to the extra dimension as 

+ E e 2ff f/gV) x \ Ln (y) + x \ R n{y)) , 

E jR (x»,y) = —L^e 2 °E { ^) X] Ro{y) 

+ -7^ E ^ 2ct xf^Cy) + £j?0O xf^y)) • (60) 

V2ttR n=0 \ J 

Here the zero mode leptons Xnoiv) and XjRoiy) are given in eq.f loTj) with the replacements 
r — )■ and r — )■ r^, respectively. 

The zero mode fermions can get mass through the Z 2 invariant left handed fermion-right 
handed fermion-Higgs interaction, tpRipL If the SM Higgs field lives on the visible brane 
as in our choice, the masses of fermions are calculated by using the integral 



rrii = 

2tt 



I rirR 

S dy \ 5 XiLo(y)XiRo(y) < H > 5(y - tt R) , (61) 

T R J-ttR 



where A5 is the coupling in five dimensions and it can be parametrized in terms of the one in four 
dimensions, the dimensionless coupling A, A5 = X/y/k. Here the expectation value of the Higgs 



10 Here, we consider different location parameters r for each left handed and right handed part of different 
flavors. The location parameters for fermion fields are chosen so that this interaction creates the current masses 
of fermions. 
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field < H > reads < H >= v/y/k where v is the vacuum expectation value n l. Now, we choose 
two different sets of location of charged lepton fields in order to obtain the masses of different 
flavors. In Set I, the left and right handed fields of the same flavor have the same location, 





SET I 


SET II 




r L r R 


r L r R 


e 


0.6710 0.6710 


-0.4900 0.8800 


n 


0.5826 0.5826 


-0.4900 0.7160 


T 


0.5273 0.5273 


-0.4900 0.6249 



Table 2: Two possible locations of charged lepton fields. Here r L and r R are left handed and 
right handed lepton field location parameters, respectively. 

however, in set II, we choose the left handed charged lepton locations the same for each flavor. 
For both cases we estimate the left-right handed charged lepton locations by respecting their 
current masses. 

The Z2 the projection condition Z^tp = —75^, used in constructing the left handed fields 
on the branes, results in that the left handed zero mode appears, the left (right) handed KK 
modes appear (disappear) on the branes, with the boundary conditions due to the Dirac mass 
term in the action eq. fl27"|) : 



d 

— m\ 

Xm n (yo) = 0, (62) 



where y = or n R. Using the Dirac equation for KK mode leptons one gets the left handed 
lepton xliniv) that lives on the visible brane as 

Xlmiv) = N Ln e°/ 2 (j^ r Je°x nL ) + c L n_ rii (e CT x ni )) , (63) 

with the constant 

J- r . T -i(x nL ) 

ei =-w^- (64) 

Here, N Ln is the normalization constant and x n L = 222 £— . The functions Jp(w) and Yp(w) 
appearing in eq. f)63p are the Bessel function of the first kind and of the second kind, respectively. 
On the other hand, the Z2 projection condition Z21P = 75^ is used in order to construct the 



n We take v = 0.043 Mpi to provide the measured gauge boson masses [JO] and choose kR = 10.83 in order 
to get the correct effective scale on the visible brane, i.e., My/ = e~^ kR M„i is of the order of TeV. 
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right handed fields on the branes and this ensures that the right handed zero mode appears, the 
right (left) handed KK modes appear (disappear) on the branes with the boundary conditions: 



f -r + m) x f Rn {y ) = o 



dy 



Again, using the Dirac equation for KK mode leptons, one gets the right handed lepton xfitniv) 
that lives on the visible brane as 

xf Rn (y) = N Rn e^ 2 U h+riH (e a x nR ) + c R Yi +nR {e° x nR )j , (66) 

with 

Cr = ~y — 77^' (67) 

where N Rn is the normalization constant and x nR = ^jr 2 -- Notice that the constant Cl, the n th 
KK mode mass m^ n in eq. (l63]) and the constant c R , the n th KK mode mass m Rn in eq. (l66|) are 
obtained by using the boundary conditions eq.( l6"2~j) and eq. (l6"5|) . respectively. For m L m\ n <C k 
and kR 3> 1 they are approximated as: 

m Ln ~ kir(n-±^ + ^)e-* kR , 

m Rn ~ kix (n - - + -) e~ nkR for r< 0.5, 

m Rn ~ kn{n+^ -) e" KkR for r> 0.5. (68) 
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^ (m) ^ (m) 



I 1 / k 



(n) 



Figure 1: One loop diagrams contribute to l\ — > / 2 7 decay due to the zero mode (KK mode) 
leptons and Higgs fields in the 2HDM. 
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Figure 2: k dependence of the BR( / u -> ei) for £f jTe = 0.01 GeV, £f iTM = 1.0 GeV. Here the 
solid (dashed, short dashed) line represents the BR without KK modes of leptons and new 
Higgs bosons (with KK modes of leptons and new Higgs bosons for lepton location set II, set 
I), for a = 0.01 and 0.1. 




Figure 3: a dependence of the BR(/i — > e^y) for £§ Te = 0.01 GeV, £f jT/i = 1.0 GeV, for the 
lepton location set II and k = 10 18 GeV. 



24 



cq 

x 




0.1 i 



0.01 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

10 18 x k (GeV) 



Figure 4: The same as Fig [3] but for r — > ej decay and for £§ Te = 0.1 GeV, £§ TT = 50 GeV. 
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Figure 5: The same as FigJ3]but for r — > ej decay, for £§ Te = 0.1 GeV, £,§ TT = 50 GeV and 
k = 10 17 GeV. 
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Figure 6: The same as Figf2]but for r — > /ry decay and for £jv T „ = 1-0 GeV, Cn,tt = §QGeV 
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Figure 7: The same as Figl3]but for r — > /I'j decay, for £jf T „ = 1.0 GeV, £§ TT = 50 GeV and 
k = 10 17 GeV. 
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